Dynamic Phase Transitions in Cell Spreading 



O 

o 



Hans-Giinthcr Dobereiner, Benjamin Dubin-Thaler, Gregory Giannone, Harry S. Xenias, and Michael P. Sheetz 
Department of Biological Sciences, Columbia University, New York, NY 10027^ 

We monitored isotropic spreading of mouse embryonic fibroblasts on fibronectin-coated substrates. 
Cell adhesion area versus time was measured via total internal reflection fluorescence microscopy. 
Spreading proceeds in well-defined phases. We found a power-law area growth with distinct expo- 
nents ai in three sequential phases, which we denote basal (ai = 0.4±0.2), continous (a-z = f .6±0.9) 
and contractile (03 = 0.3 ± 0.2) spreading. High resolution differential interference contrast mi- 
croscopy was used to characterize local membrane dynamics at the spreading front. Fourier power 
spectra of membrane velocity reveal the sudden development of periodic membrane retractions at 
the transition from continous to contractile spreading. We propose that the classification of cell 
spreading into phases with distinct functional characteristics and protein activity patterns serves as 
a paradigm for a general program of a phase classification of cellular phenotype. Biological vari- 
ability is drastically reduced when only the corresponding phases are used for comparison across 
species/different cell lines. 
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PACS numbers: 05.45.-a, 87.17.-d 



Cells need to be mobile in order to perform many crit- 
ical biological functions. The reorganization of extracel- 
lular matrix in wound healing, the positioning of nerve 
cells, or the engulfment of bacteria in the immune reac- 
tion of white blood cells are particular examples Q] . Ac- 
complishing this variety of functions requires a diverse 
set of mechanisms and proteins. Most components of the 
molecular machinery of actin-based motility have been 
identified 0,111 0]- It has been possible to perform exper- 
iments with reconstituted systems of Listeria propulsion 
[jj |(| for which detailed elastic models have been devel- 
oped 0. In vivo, cell spreading on matrix-coated surfaces 
provides a simplified system of analyzing motile behav- 
ior. A substantial amount of experimental and theoret- 
ical work has been done along these lines H 0, 0, Q] . 
However, only quite recently, quantitative experiments 
of whole cell spreading and subsequent migration could 
be performed with high spatial and temporal resolution 
[l2l If3| . We found that there are well defined and dis- 
tinct states of spreading. It is the goal of this work to 
show that these states can indeed be considered phases 
of motility by demonstrating the existence of dynamic 
phase transitions between them. 

Spreading cells extend a 200 nm thick sheet called the 
lamellipodium from the cell body onto the substrate, see 
Fig. ^ This process is driven by actin polymerization 
at the leading membrane edge, the precise mechanism 
of which is still under debate [fj, [l5j. The meshwork 
of actin fibers is crosslinked by various proteins. The 
molecular motor myosin II enables the meshwork to con- 
tract by moving along actin fibers and relative to other 
cytoskclctal elements. Thus, the lammelipodium is an 
active gel enclosed in a flat membrane bag adhering to 
the substrate. The physics of active gels has recently 
attracted a lot of attention. Rheological experiments of 
simple mixtures of purified actin and myosin solutions 
[lfij l and quite general theoretical modeling 0, 0] have 



been carried out. There are dynamic phase transitions 
involving extended and contracted actin density states as 
a function of myosin-actin coupling strength |l7| . We will 
show that our cellular system exhibits similar transitions 
which express themselves prominently in the dynamics of 
the leading membrane edge. 




FIG. 1: A: During cell spreading, a thin lamellipodial sheet 
extends from the cell body onto the substrate. B: Total in- 
ternal reflection fluorescence micrograph of a spreading cell. 
The bright region corresponds to the area adhered to the sub- 
strate. C: Two overlayed snapshots of the leading membrane 
edge of a lamellipodium moving from right to left are shown in 
differential interference contrast. The edge position is marked 
with a white contour overlay. 



Mouse embryonic fibroblasts (MEF) were allowed to 
settle onto fibronectin coated glass slides and observed 
with either total internal reflection fluorescence (TIRF) 
or differential interference contrast (DIC) microscopy. 
Fibronectin links the cell membrane to the extracellular 
matrix proteins mimicking natural cell adhesion. TIRF 
studies were performed at a moderate spatial and tempo- 
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ral resolution to capture overall spreading characteristics 
of the whole cell. Multiple cells could be studied simul- 
taneously. High resolution DIC was used to characterize 
local membrane dynamics. Details of the methods may 
be found in our earlier work |l2l . 
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FIG. 2: Adhesion area in isotropically spreading fibroblasts 
grows with a power law in time. Different but constant ex- 
ponents Oi in the various phases of spreading are evident in 
a double logarithmic plot. Exponents have been determined 
by fitting a piecewise linear function to the data, see Fig. [3] 
Adhesion areas Ai at the transition points are indicated. 



Membrane adhesion area A during spreading was 
best monitored using T1RF. Distinct classes of angular 
isotropic and anisotropic spreading cells were found |l2| . 
In the following, we limit ourselves to the isotropic class 
lacking filopodia. Close inspection of double logarithmic 
plots of adhesion area A over time reveals three phases 
with distinctly different power law growth, as seen in 
Fig. |2] We define area growth exponents a* via 



A(t) ~ i* , 



(1) 



where i denotes the subsequent phases. Initially, there is 
a basal phase where cells test the suitability of the sub- 
strate to adhere and area growth is mimimal. We find 
ai = 0.4 ± 0.2. Then follows a phase of fast continuous 
spreading, which is characterized by a 2 = 1.6 ± 0.9. Fi- 
nally, the cell slows down again exhibiting a sub-linear 
area growth with 0,3 = 0.3 ± 0.2. We will see below that 
the latter phase is characterized by periodic local con- 
tractions of the cell Nevertheless, the mean area 
growth leads to an effective power law behavior also in 
this phase. Histograms of exponents Oj for the three 
phases are shown in Fig. [3J There is a clear distinction of 
fast area growth in the continuous spreading phase with 
a rather broad distribution of the exponent a-i- How- 
ever, we find two narrow clusters when discriminating 
with respect to the relative area growth, A2/A1, during 



that phase, where A4 denotes the adhesion area at the 
transition from phase i to i + 1. Small (A2/A1 < 5) 
or large (A2/A1 > 5) area increases correspond to small 
(rn2 = 0.9 ± 0.2) or large (77J2 = 1.6 ± 0.2) exponents, 
respectively. In addition, there were two single cells with 
even larger exponents mi which we excluded from the 
cluster average. 
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FIG. 3: Histograms of area growth exponents as obtained 
from the slopes of double logarithmic plots of adhesion area 
versus time, like the ones shown in Fig. |21 We have analyzed 
20 cells in total. The middle (continous spreading) phase 
exhibits clustering corresponding to small (open bars) and 
large (filled bars) area growth during that phase. 



The transition from continuous to contractile spread- 
ing was further monitored using high resolution DIC. 
A suitable isotropically spreading cell was chosen and 
a well- resolvable and approximately straight membrane 
segment was selected for prolonged observation, see 
Figffl B. Time-lapse sequences were obtained at video 
rate. Movies were digitized at 1/At = 3Hz. Individual 
frames are counted using an index n. The cell edge is 
determined with a custom C program by a local contour 
algorithm 19] allowing nanometer accuracy. We obtain a 
sub-pixel resolution of 15 nm for relative displacements, 
which translates into a minimal detectable velocity of 45 
nm/s between frames. Further analysis proceeds using 
a cartesian coordinate system where the average mem- 
brane orientation is taken as the fixed y-axis. Points on 
the membrane are then labeled by their y-coordinates yj 
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and the membrane velocity Vj(n) = Axj(n)/ At is mea- 
sured along the x-axis which is normal to the average 
membrane orientation. 
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FIG. 4: Normal velocity map of the particular membrane 
segment marked in Fig. 03 as a function of time. Note the two 
qualitatively different sections before and after time t = 40 s 
corresponding to a continous and a periodically contractile 
spreading phase, respectively. The period of the latter is T = 
17±4 s. The speeds of lateral waves of maximum contraction 
velocity are indicated. The encircled region marks a phase 
shift of a contraction. 



A typical velocity map along the contour over time is 
shown in Fig. 0] We find that a region of continous, un- 
interrupted spreading (red shadows) precedes a sequence 
of periodic membrane retraction events (blue stripes). 
These two different states of membrane dynamics corre- 
spond to the continuous spreading and contractile phase 
of the lamcllipodium, found above. The two phases can 
be clearly distinguished using the discrete Fourier trans- 
formation v(s) of the velocity map v(n) defined as 

1 / >. /„ (n-l)(s-l)\ 
v i ( s)= N^ v i ( n ) ex P 2m ~ S " ) > ( 2 ) 

71=1 ^ ' 

where N is the number of frames. Averages are taken 
over spatial regions of interest. The continuous spread- 
ing phase is characterized by a strong boundary maxi- 
mum of the spectrum |vj(s)| at s = 1, see Fig. 0A. In 
contrast, in the contractile phase the spectrum develops a 
pronounced peak at s = s max , see Fig.[S]B, which signals 
oscillatory behavior with a period 

N 

T = At . (3) 

Thus, the peak position of the power spectrum serves 
as an excellent phase indicator. We calculate the spec- 
trum inside a small time window - with a width on the 
order of the repeat time - and sweep across the phase 
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FIG. 5: Fourier spectrum (see Eq. of the velocity map in 
Fig. 0| for the two different spreading phases below (A,N = 
120) and above (B,N = 200) frame number n = 120. The 
spectrum is spatially averaged over 70 points between posi- 
tion 2.0 fim and 3.8 (im along the contour. The transition 
between continous and contractile spreading is characterized 
by a sharp shift in the position of the maximum of the Fourier 
spectrum. The boundary maximum in panel A corresponds to 
a mean velocity of 7.4 /im in the continous spreading phase. 
The peak at s max = 5 in panel B corresponds to a repeat 
time T — 17 ± 4 s for the contractile spreading phase, see 
Eq.|H] The mean velocity 0.3 /im is small. The lower panel C 
depicts the peak position of the spectrum taken in a running 
time window with a width of N — 50 frames, corresponding 
to the repeat time T, as a function of the first frame number 
of the window. Note that the peak position depends on the 
width of the window used for Fourier transformation. 



boundary. Indeed, there is a well defined transition be- 
tween the two phases as seen in Fig. 0C. However, the 
periodic contractions do not take place simultaneously 
along the leading edge, see Fig.@] In fact, there are lat- 
eral waves of maximum contraction velocity running in 
both directions. These waves have a speed on the order 
of 200 nm/s. Moreover, there are sudden phase shifts of 
the periodic contractions up to half a period, see encir- 
cled region in Fig. 0] 

To summarize our experimental findings, we have seen 
clear signatures of dynamic phase transitions in the 
spreading behavior of MEF cells. Since actin poly- 
merization does not stop during membrane retraction 
events ^3[> one concludes that the actin network con- 
tracts and/or is actively pulled back by myosin motor 
activity. Kruse et al. [lj| have modeled oriented fibers 
connected into a network by molecular motors. They 
find an instability of homogeneous fiber density towards 
a contracted state as a function of fiber-motor coupling 
strength. Moreover, their generic theoretical model al- 
lows for oscillatory solutions. Our cellular system ex- 
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hibits similar behavior. Indeed, the periodic contrac- 
tions are absent when myosin light chain kinase (MLCK) , 
which activates myosin, is inhibited |l3| . 

In the following, we give a systems biology oriented 
view of cell motility. Several questions arise: What is the 
functional role of the phases described above? How are 
these dynamic phases of the structural motility network 
regulated and connected to the signaling network? Can 
we disentangle the complex set of motility related pro- 
teins and simplify description by considering functional 
modules and hierarchical levels of control? 

During the initial spreading phase, a MEF cell assem- 
bles the cytoskeletal structure necessary to probe the 
mechanical suitability of the substrate in the following 
contractile phase where it periodically pulls on the sub- 
strate. Indeed, cells require stiff substrates for growth 
and move toward stiffer regions fll| . The machinery for 
this stiffness sensing is organized into i) the basic struc- 
tural elements consisting of the actin cytoskeleton, the 
myosin II motors, as well as the plasma membrane, ii) 
factors directly regulating protein coupling strength and 
activity, and, finally, iii) these regulatory proteins are 
controlled by a signaling network coordinating spatially 
distant and/or logically separate functional events in the 
cell. In order to link these cellular components to the 
physics of dynamic phase transitions, we note that it is 
the basic structural elements which exhibit the various 
phases we have identified in this work. The phase param- 
eters are given by the regulating proteins, e.g., MLCK ac- 
tivity and concentration. The different dynamic phases 
correspond to functional regions in the regulating para- 
meter space. The trajectories in this parameter space are 
determined by the cellular signaling network. This hier- 
archical identification provides an immediate conceptual 
advantage: The topology of the motility phase diagram 
is independent of the complex signaling network, i.e., the 
relative positions of all the motility phases do not de- 
pend on the trajectories in parameter space followed by 
the cell. Indeed, the basic phase characteristics can be 
probed and modeled separately. This was demonstrated 
using the observed linear relationship between the period 
of contractions and the lamellipodia width 0] . We find 
that the contraction period is the same for equal lamel- 
lipodia width, independent from the variations in the bio- 
chemical pathway(s) induced in order to achieve a certain 
width. On the other hand, cell motility is a unique case 
where the interconnections of the structural and signal- 
ing networks, which depend on gene expression, can be 
probed in order to establish a quantitative link between 
phenotype and genotype. We are currently identifying 
functional modules in a large scale screening of spread- 
ing phenotypes across various fibroblast cell lines with 
mutant genotypes. 

The idea of phases in cell behavior can be applied 
quite generally. Phases of motility should be consid- 
ered analogous to the phases of the cell cycle, phases of 



varying metabolic activity or different protein expression. 
We propose to classify cellular behavior in well defined 
phases. Their number will be considerable less than an 
enumeration of concentration and activity levels of all 
molecular components of the cell. Thus, one can hope 
to accomplish a simplified description. Currently, phase 
classification is not generally done and cellular phenotype 
cannot be sensibly compared across different genotypes. 
We expect that some fraction of the variability encoun- 
tered in biological experiments and the often conflicting 
results between laboratories stem from the fact that find- 
ings corresponding to different cellular phases and bound- 
ary conditions are spuriously compared to each other. In 
conclusion, we feel that the classification of motility in 
phases can serve as a paradigmatic example for a power- 
ful general ordering principle in quantitative biology. 
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